clear all

use "Estimates_for_plots_6_new_reference.dta"

destring month, replace
rename estimate coef


graph set window fontface "Times New Roman"


**AUSTRIA, 2009
*++++++++++++++++++++
**Subsidized and non-subsidized unit sales, time path, window around reform [-6,6]
twoway	(rarea max95 min95 month if type=="s" & program=="at09", color(navy*.3)) (rarea max95 min95 month if type=="ns" & program=="at09", color(none) lcolor(black) lpattern(dash) msize(small))(scatter coef month if type=="s" & program=="at09", connect(l) lpattern(solid) lcolor(black) msize(small) mcolor(black) msym(D) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) graphregion(color(white))  legend(off) xtitle("Month")  ytitle("log(UNITS)") ylabel(-1.0(.5)1.0)) (scatter coef month if type=="ns" & program=="at09" , connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Mar." -5 "Apr." -4 "May" -3 "Jun."-2 "Jul." -1 "Aug." 0 "Sept." 1 "Oct." 2 "Nov." 3 "Dec." 4 "Jan." 5 "Feb." 6 "Mar." 7 "Apr."8 "May" 9 "Jun." )) 

graph export "AT_09_u_62.eps", as(eps) fontface("Times New Roman") replace


**Subsidized and non-subsidized prices, time path, window around reform [-6,6]
twoway	(rarea max95 min95 month if type=="sp" & program=="at09" , color(navy*.3)) (rarea max95 min95 month if type=="nsp" & program=="at09", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="sp" & program=="at09" , connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1)) (scatter coef month if type=="nsp" & program=="at09", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Mar." -5 "Apr." -4 "May" -3 "Jun."-2 "Jul." -1 "Aug." 0 "Sept." 1 "Oct." 2 "Nov." 3 "Dec." 4 "Jan." 5 "Feb." 6 "Mar." 7 "Apr."8 "May" 9 "Jun." ))

graph export "AT_09_p_62.eps", as(eps) fontface("Times New Roman") replace




**AUSTRIA, 201O, FF
*++++++++++++++++++++
**Subsidized and non-subsidized unit sales, time path, window around reform [-6,6]
twoway(rarea  max95 min95 month if type=="s" & program=="at10ff" , color(navy*.3)) (rarea max95 min95 month if type=="ns" & program=="at10ff", color(none) lcolor(black) lpattern(dash))(scatter coef month if type=="s" & program=="at10ff", connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1) ) tline(2, lstyle(refline) lwidth(0.1) lcolor(gs8)  lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1) ) lpattern(solid)lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(UNITS)") ylabel(-1.0(.5)1.0)) (scatter coef month if type=="ns" & program=="at10ff", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Mar." -5 "Apr." -4 "May" -3 "Jun."-2 "Jul." -1 "Aug." 0"Sept." 1 "Oct." 2 "Nov." 3"Dec." 4 "Jan." 5 "Feb." 6 "Mar." 7 "Apr."8 "May"))

graph export "AT_10_ff_u_62.eps", as(eps) fontface("Times New Roman") replace

**Subsidized and non-subsidized prices, time path, window around reform [-6,6]
twoway(rarea max95 min95 month if type=="sp" & program=="at10ff" , color(navy*.3)) (rarea max95 min95 month if type=="nsp" & program=="at10ff", color(none)lcolor(black)lpattern(dash))(scatter coef month if type=="sp" & program=="at10ff", connect(l)tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1) ) tline(2, lstyle(refline) lwidth(0.1) lcolor(gs8)  lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1) ) lpattern(solid)lcolor(black) mcolor(black) msym(D)msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(PPRICE)") ylabel(-0.05(.05)0.1) )(scatter coef month if type=="nsp" & program=="at10ff", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Mar." -5 "Apr." -4 "May" -3 "Jun."-2 "Jul." -1 "Aug." 0"Sept." 1 "Oct." 2 "Nov." 3"Dec." 4 "Jan." 5 "Feb." 6 "Mar." 7 "Apr."8 "May"))

graph export "AT_10_ff_p_62.eps.eps", as(eps) fontface("Times New Roman") replace



**HUNGARY, 2016
*++++++++++++++++++++
****Subsidized and non-subsidized unit sales, time path, window around reform [-6,1]. Longer window impossible, given data.
twoway(rarea  max95 min95 month if type=="s" & program=="hu16" & label=="all" , color(navy*.3)) (rarea max95 min95 month if type=="ns" & program=="hu16" & label=="all", color(none) lcolor(black) lpattern(dash))(scatter coef month if type=="s" & program=="hu16" & label=="all", connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid)lcolor(black) mcolor(black) msym(D)msize(small) graphregion(color(white)) legend(lab(3 "Time Path of Unit Sales Subsidized") lab(4 "Time Path of Unit Sales Non-Subsidized") lab(1 "95% CI") lab(2 "95% CI") region(lstyle(none) color(none)) nobox) leg(on order(3 1 4 2 ))xtitle("Month") ytitle("log(UNITS)") ylabel(-1(.5)1))(scatter coef month if type=="ns" & program=="hu16" & label=="all", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Mar." -5 "Apr." -4 "May" -3 "Jun."-2 "Jul." -1 "Aug." 0"Sept." 1 "Oct." 2 "Nov." 3"Dec." 4 "Jan."))

graph export "HU_16_u_62.eps", as(eps) fontface("Times New Roman") replace

****Subsidized and non-subsidized prices, time path, window around reform [-6,1]. Longer window impossible, given data.
twoway(rarea max95 min95 month if type=="sp" & program=="hu16"& label=="all" , color(navy*.3)) (rarea max95 min95 month if type=="nsp" & program=="hu16"& label=="all", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="sp" & program=="hu16"& label=="all" , connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(gs8)  lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(lab(3 "Time Path of Price Subsidized") lab(4 "Time Path of Price Non-Subsidized") lab(1 "95% CI") lab(2 "95% CI") region(lstyle(none)color(none)) nobox) leg(on order(3 1 4 2 )) xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1)) (scatter coef month if type=="nsp" & program=="hu16" & label=="all", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Mar." -5 "Apr." -4 "May" -3 "Jun."-2 "Jul." -1 "Aug." 0"Sept." 1 "Oct." 2 "Nov." 3"Dec." 4 "Jan." ))

graph export "HU_16_p_62.eps", as(eps) fontface("Times New Roman") replace


**AUSTRIA, 2010, WM
*++++++++++++++++++++
**Subsidized and non-subsidized unit sales, time path, window around reform [-6,6]
twoway(rarea max95 min95 month if type=="s" & program=="at10wm" , color(navy*.3)) (rarea max95 min95 month if type=="ns" & program=="at10wm", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="s" & program=="at10wm", connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(1, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(UNITS)") ylabel(-1(.5)1.0)) (scatter coef month if type=="ns" & program=="at10wm", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Oct." -5 "Nov." -4 "Dec." -3 "Jan."-2 "Feb." -1 "Mar." 0"Apr." 1 "May." 2 "Jun." 3"Jul." 4 "Aug." 5 "Sep." 6 "Oct." 7 "Nov."))

graph export "AT_10_wm_u_62.eps", as(eps) fontface("Times New Roman") replace

**Subsidized and non-subsidized prices, time path, window around reform [-6,6]
twoway(rarea max95 min95 month if type=="sp" & program=="at10wm" , color(navy*.3)) (rarea max95 min95 month if type=="nsp" & program=="at10wm", color(none) lcolor(black) lpattern(dash))(scatter coef month if type=="sp" & program=="at10wm", connect(l)tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(1, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1)) (scatter coef month if type=="nsp" & program=="at10wm", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Oct." -5 "Nov." -4 "Dec." -3 "Jan."-2 "Feb." -1 "Mar." 0"Apr." 1 "May." 2 "Jun." 3"Jul." 4 "Aug." 5 "Sep." 6 "Oct." 7 "Nov."))

graph export "AT_10_wm_p_62.eps", as(eps) fontface("Times New Roman") replace



**HUNGARY, 2015
*++++++++++++++++++++
**Subsidized and non-subsidized unit sales, time path, window around reform [-6,6]
twoway(rarea max95 min95 month if type=="s" & program=="hu15" & label=="all" , color(navy*.3)) (rarea max95 min95 month if type=="ns" & program=="hu15" & label=="all", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="s" & program=="hu15" & label=="all", connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(lab(3 "Time Path of Unit Sales Subsidized") lab(4 "Time Path of Unit Sales Non-Subsidized") lab(1 "95% CI") lab(2 "95% CI") region(lstyle(none) color(none)) nobox) leg(on order(3 1 4 2 )) xtitle("Month") ytitle("log(UNITS)") ylabel(-1(.5)1))(scatter coef month if type=="ns" & program=="hu15" & label=="all", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Apr." -5 "May" -4 "Jun."-3 "Jul."-2 "Aug." -1 "Sep." 0"Oct." 1 "Nov." 2 "Dec." 3"Jan." 4 "Feb." 5 "Mar." 6 "Apr." 7 "May" 8 "Jun." 9 "Jul." ))

graph export "HU_15_u_62.eps", as(eps) fontface("Times New Roman") replace

**Subsidized and non-subsidized prices, time path, window around reform [-6,6]
twoway(rarea max95 min95 month if type=="sp" & program=="hu15" & label=="all" , color(navy*.3)) (rarea max95 min95 month if type=="nsp" & program=="hu15" & label=="all", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="sp" & program=="hu15" & label=="all", connect(l)tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1) ) tline(3, lstyle(refline) lwidth(0.1) lcolor(gs8) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D)msize(small) graphregion(color(white)) legend(lab(3 "Time Path of Price Subsidized") lab(4 "Time Path of Price Non-Subsidized") lab(1 "95% CI") lab(2 "95% CI") region(lstyle(none) color(none)) nobox) leg(on order(3 1 4 2 )) xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1)) (scatter coef month if type=="nsp" & program=="hu15" & label=="all", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black)  msym(O) xlabel(-6 "Apr." -5 "May" -4 "Jun." -3 "Jul."-2 "Aug." -1 "Sep." 0"Oct." 1 "Nov." 2 "Dec." 3"Jan." 4 "Feb." 5 "Mar." 6 "Apr." 7 "May" 8 "Jun." 9 "Jul."))

graph export "HU_15_p_62.eps", as(eps) fontface("Times New Roman") replace



**CROATIA, 2015, WM
*++++++++++++++++++++
**Subsidized and non-subsidized unit sales, time path, window around reform 
twoway(rarea max95 min95 month if type=="s" & program=="hr15" , color(navy*.3)) (rarea max95 min95 month if type=="ns" & program=="hr15", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="s" & program=="hr15", connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(4, lstyle(refline) lwidth(0.1) lcolor(gs8)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(UNITS)") ylabel(-1(.5)1)) (scatter coef month if type=="ns" & program=="hr15", connect(l) lpattern(dash) lcolor(black) msize(small) mcolor(black) msym(O) xlabel(-6 "Dec." -5 "Jan." -4 "Feb." -3 "Mar." -2 "Apr." -1 "May" 0"Jun." 1"Jul." 2 "Aug." 3 "Sept." 4"Oct." 5"Nov." 6 "Dec." 7 "Jan." 8 "Feb." 9 "Mar." 10 "Apr."))

graph export "HR_15_u_62.eps", as(eps) fontface("Times New Roman") replace

**Subsidized and non-subsidized prices, time path, window around reform 
twoway(rarea max95 min95 month if type=="sp" & program=="hr15" , color(navy*.3)) (rarea max95 min95 month if type=="nsp" & program=="hr15", color(none) lcolor(black) lpattern(dash)) (scatter coef month if type=="sp" & program=="hr15", connect(l) tline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) tline(4, lstyle(refline) lwidth(0.1) lcolor(gs8)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(off) xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1)) (scatter coef month if type=="nsp" & program=="hr15", connect(l) lpattern(dash) lcolor(black) msize(small)  mcolor(black) msym(O) xlabel(-6 "Dec." -5 "Jan." -4 "Feb." -3 "Mar." -2 "Apr." -1 "May" 0"Jun." 1"Jul." 2 "Aug." 3 "Sept." 4"Oct." 5"Nov." 6 "Dec." 7 "Jan." 8 "Feb." 9 "Mar." 10 "Apr."))

graph export "HR_15_p_62.eps", as(eps) fontface("Times New Roman") replace

